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This  report  describes  the  Naval  Postgraduate  School  (NPS)  Environmental 
Physics  Group's  approach  to  predicting  properties  of  the  marine  atmospheric 
boundary  layer  which  affect  weapons  systems.  We  describe  tactically  signi¬ 
ficant  parameters  and  how  they  can  be  predicted  using  a  boundary  layer  model. 
We  also  describe  verification  of  the  model  with  data  gathered  from  shipboard 
sensors.  This  is  the  second  in  a  series  of  reports  on  the  NPS  marine  atmo¬ 
spheric  prediction  model;  the  first  report  was  a  position  paper  on  the  model¬ 
ing  method  (Fairall  et  al ,  1981). 

Modern  warfare  has  become  critically  dependent  on  nearly  the  entire 
electromagnetic  spectrum  for  command,  control  and  communications  for  weapon 
guidance,  for  electronic  warfare  support  and  for  countermeasures.  Thus,  en¬ 
vironmental  enhancement  or  degradation  of  the  performance  of  electromagnetic 
(EM)  and  electro-optic  (EO)  systems  has  become  a  primary  concern  of  task  force 
commanders.  Tactically  essential  systems  are  affected  by  the  environment  even 
when  meteorological  conditions  are  not  severe  in  the  historical  sense.  The 
deployment  of  systems  and  the  modification  of  tactics  based  on  environmental 
factors  will,  to  a  large  extent,  determine  the  effectiveness  of  EM/EO  sensor, 
weapon  and  communication  systems. 

The  assessment  of  environmental  effects  on  EM/EO  systems  is  not  a  simple 
task  because  there  are  many  variables  to  consider.  A  given  platform  may  have 
many  systems,  each  with  its  individual  characteristics  and  purpose,  and  per¬ 
formance  depends  on  such  systems  parameters  as  frequency,  power,  antenna  gain 
and  on  such  meteorological  and  oceanic  parameters  as  wind  speed,  humidity, 
temperature  and  sea  state.  If  a  complete  and  accurate  set  of  parameters  de¬ 
scribing  the  environment  is  available,  the  environmental  impact  on  weapon 
systems  can  usually  be  calculated.  However,  the  task  force  commander  often 
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has  an  incomplete  set  of  environmental  parameters  available  to  him  and  in  a 
wartime  situation  his  communication  with  external  sources  of  information,  such 
as  Fleet  Numerical  Oceanography  Center,  may  even  be  cut.  It  is  therefore 
necessary  for  task  force,  or  even  individual  ship  commanders  to  be  able  to 
make  spatial  and  temporal  extrapolations  of  environmental  parameters  based 
solely  on  information  they  themselves  can  obtain  and  computational  capabili¬ 
ties  that  reside  with  them. 

A  gap  has  existed  in  past  efforts  to  characterize  tactical  environmental 
conditions.  The  gap  was  between  the  two  extreme  approaches  of  relating  condi¬ 
tions  to  1)  near  surface  observations,  and  to  2)  larger  scale  predicted  synop¬ 
tic  patterns.  Clearly,  to  assess  the  local  features,  local  measurement  is 
desirable.  However,  to  assess  over  a  large  spatial  region  or  for  a  long  time 
period,  a  local  measurement  made  at  one  time  loses  its  validity.  Thus,  one 
must  consider  a  transition  to  climatology,  or  large  scale  numerical  analysis 
predictions,  or  dynamic  model s  based  on  observations  at  the  operational  loca¬ 
tion.  Climatologies  do  not  include  the  small  scale  structure  descriptions 
needed.  Ruggles  (1975)  has  examined  the  capabilities  of  large  scale  numerical 
procedures  at  remote  facilities  and  argued  convincingly  that  they  are  not 
sufficient  for  all  tactical  descriptions.  Thus,  we  must  rely  on  dynamic  models 
that  utilize  both  local  measurements  and  large  scale  numerical  predictions. 
This  is  the  philosophy  behind  the  Tactical  Environmental  Support  System  (TESS) 
concept  and  the  NPS  model  we  describe  here. 

Three  types  of  models  are  needed  to  predict  weapon  system  behavior.  One 
is  needed  to  calculate  system  performance  given  a  set  of  environmental  para¬ 
meters.  An  example  is  the  Integrated  Refractive  Effects  Prediction  System 
(IREPS)  which  determines  and  displays  parameters  affecting  radar  and  communi¬ 
cation  systems  using  mean  meteorological  parameters  as  input  information. 
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Another  is  required  to  calculate  environmental  parameters  from  readily  avail¬ 
able  meteorological  data.  An  example  is  the  estimation  of  thermal  turbulence 
needed  to  assess  optical  system  performance,  which  cannot  be  measured  directly 
but  can  be  determined  from  radiosonde  data.  Finally,  another  is  needed  to 
predict  the  changes  in  the  meteorological  properties. 

It  is  assumed  at  the  outset  that  models  should  be  developed  which  depend 
as  much  as  possible  on  "single  station  assessments".  That  is,  they  should 
only  require  data  which  can  be  gathered  from  a  single  ship  or  ation.  The 
models  should  also  have  the  capability  to  use  data  from  many  trees,  which 
would  improve  the  accuracy  of  the  assessment,  but  multi-source  a  would  not 
be  a  necessity.  Another  requirement  is  that  it  should  be  possi,  •-  o  run  the 
model  on  a  microcomputer,  which  would  make  it  available  to  a  wide  range  of 
users  and  compatible  with  the  microcomputer  based  IREPS  program. 
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2.  ATMOSPHERIC  MODELING  FOR  TACTICAL  DESCRIPTIONS 

Predicting  the  effect  of  the  environment  on  tactical  systems  ranges  from 
determining  large  scale  atmospheric  forcing  on  the  boundary  layer  to  assess¬ 
ment  of  the  behavior  of  an  actual  system.  The  latter  requires  the  placement 
of  a  set  of  instructions  in  the  hands  of  the  systems  operator.  This  section 
describes  the  requirements  on  a  boundary  layer  model,  which  is  to  be  used  to 
predict  changes  in  significant  propagation  parameters  within  the  first  kilo¬ 
meter.  The  various  models  and  their  purposes  within  the  tactical  situation 
are  tested  in  Table  1.  The  boundary  layer  model  is  only  one  of  the  models;  it 
calculates  the  expected  evolution  of  the  boundary  layer  properties  and  outputs 
this  information  for  use  by  propagation  models. 

A  simplified  schematic  of  the  information  flow  between  the  models  for  en¬ 
vironmental  assessment  is  shown  in  Figure  1.  The  predicted  properties  would 
be  used  by  such  assessment  systems  as  IREPS  for  radar  and  radio  propagation, 
LOWTRAN  and  SAEL  for  optical  propagation.  The  latter  are  shown  on  the  RHS  of 
Figure  1. 

The  two  primary  properties  in  the  assessment  of  electromagnetic  propaga¬ 
tion  are  the  index  of  refraction  and  extinction  parameters. 

Extinction:  Extinction  depends  on  the  electromagnetic  wavelength.  In  the 
radio  and  radar  region  extinction  is  due  to  scattering  by  droplets  and  from 
the  sea  surface  and,  to  a  lesser  extent,  scattering  from  refractive  index  in¬ 
homogeneities.  In  general,  scattering  is  a  minor  effect  on  radio  and  radar 
frequencies  for  most  cases  except  rain  and  high  sea  state  surface.  In  the 
optical  region,  extinction  is  due  to  molecular  absorption  and  aerosol  scatter¬ 
ing  and  can  be  quite  severe.  Molecular  absorption  occurs  over  narrow  wave¬ 
length  bands  so  that  optical  "windows"  exist.  Extinction  in  fog  and  clouds  is 
so  severe  that  most  optical  systems  cannot  be  used. 


Refraction:  Index  of  refraction  gradients  exist  over  scales  from  kilome¬ 


ters  to  centimeters.  Gradients  over  the  large  scale  lead  to  ray  path  curva¬ 
ture  which  seriously  affects  the  performance  of  radio  and  radar  systems,  but 
not  optical  systems  which  normally  are  operated  over  short  paths.  However, 
small  scale  index  variations  (1  cm  to  10  m)  can  cause  significant  degradation 
in  optical  images  and  the  energy  density  that  can  be  obtained  in  an  optical 
beam. 

The  tactical  importance  of  variations  in  index  of  refraction  gradients  and 

variations  is  illustrated  in  Figures  2  and  3.  Figure  2  illustrates  refractive 

effects  on  the  radar  propagation  index  when  the  ray  curvature  is  greater  than 

the  curvature  of  the  earth.  Figure  2a  is  for  a  surface  based  radar  and  2b  is 

for  an  airborne  radar.  Degradation  of  an  optical  image  with  increasing  small 

scale  index  of  refraction  turbulence,  in  terms  of  the  index  of  refraction 

2 

structure  function  Cn  ,  is  shown  in  Figure  3.  Values  range  from  no  turbulence. 
Figure  3a,  to  moderate  turbulence,  Figure  3c.  As  the  figures  indicate,  atmo¬ 
spheric  effects  are  of  serious  concern  if  effects  such  as  those  shown  occur  a 

significant  fraction  of  the  time. 

2 

Variation  of  Cn  with  height  is  shown  in  Figure  4.  These  data  were  ob¬ 
tained  in  an  NPS  experiment  and  are  typical  of  what  is  observed  in  and  at  the 
top  of  a  well  mixed  boundary  layer.  Cp  is  large  near  the  surface  and  inver¬ 
sion.  Reference  to  Figure  2  shows  that  values  observed  at  the  inversion  lead 
to  substantial  degradation  of  optical  images. 

2 

Beam  spread  with  distance  is  shown  in  Figure  5  for  Cn  ranging  from  low  to 
moderate  values.  For  the  lower  turbulence  value  the  beam  radius  increases  by 
a  factor  of  4  from  1  km  to  4  km,  which  corresponds  to  an  order  of  magnitude 
decrease  in  energy  density. 
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TABLE  1.  Prediction/assessments  and  models  required. 


Purpose 

Required 

Synoptic  Scale 

Forcing 

Winds 

Subsidence 

Advectic 

Synoptic  Scale 
Numerical 

Local  Meteorological 
Parameters 

Inversion 

Temperature 

Humidity 

Clouds 

Winds 

Boundary  Layer 

Propagation 

Parameters 

Index  of  Refraction 
Aerosol 

Absorbing  Molecules 

Boundary  Layer 

EM/EO  Behavior 

Ray  Curvature 

Vi  Sibil i ty 

Image  Resolution 

System  Dependent 
Propagation 

Tactical  Decision 

Aids 

System  Utilization 

FI ight  Pattern 

EM/EO  Order  of  Battle 

Specialized 
Tactical  Aids 
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Processes/Properties  Environmental  Predicted  Propagation 

Models  Properties  Models 

(TESS) 


The  predicted  meteorological  properties  that  are  necessary  to  the  above 
features  and  which  are  available  from  the  NPS  boundary  layer  model  appear  in 


the  second  column  from  the  RHS  in  Figure  1  and  are: 

Inversion  height  and  jump  magnitudes: 

The  height  of  the  inversion  defines  the  depth  of  the  marine  boundary  layer. 
This  is  also  the  elevation  where  electromagnetic  trapping  occurs.  The  in¬ 
dex  of  refraction  gradient  governs  the  path  curvature  and  the  gradient  de¬ 
pends  on  the  changes  (jumps)  in  temperature  and  water  vapor  in  the  inver¬ 
sion.  The  model  outputs  the  index  of  refraction  profile  (in  M  units)  for 
use  in  the  IREPS  program.  The  boundary  layer  depth  also  defines  the  region 
through  which  materials  released  within  the  layer  are  dispersed,  a  quanti¬ 
ty  which  is  needed  for  dispersion  models.  For  example,  aerosols  generated 
at  the  sea  surface  are  trapped  within  this  layer,  which  impacts  on  optical 
extinction. 

Boundary  layer  well  mixed  properties: 

The  equivalent  potential  temperature  and  total  water  mixing  ratio  (mea¬ 
sures  of  temperature  and  moisture)  are  basic  boundary  layer  properties 
that  are  used  for  all  calculations  within  the  model.  They  are  also  needed 
for  all  of  the  propagation  models  and  are  available  as  outputs. 

Aerosol  loading,  extinction: 

The  aerosol  equilibrium  size  distribution  as  a  function  of  height  is  cal¬ 
culated  on  the  basis  of  the  ambient  relative  humidity  and  wind  speed. 
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(a)  0  (no  turbulence) ,  ( 
lO-14  nf'^  (from  Keams 
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Figure  5.  Variation  of  laser  beam  radius  with  distance  for  Cn^  levels 
of  6. 5  x  ID'15  m  (solid  line)  and  6.5  x  10  ^  m 


A 


(dashed  line). 


Fluxes: 


The  fluxes  of  momentum,  heat  and  water  vapor  and  the  attendant  turbulence 
intensities  are  all  possible  output  parameters.  The  surface  layer  fluxes 
are  needed  to  calculate  evaporation  duct  properties.  Turbulence  intensi¬ 
ties  through  the  boundary  layer  are  needed  to  calculate  optical  scintil¬ 
lation,  wave  front  distortion  and  slant  path  propagation  properties. 

Fog  and  cloud  prediction: 

An  integral  part  of  the  model  is  the  calculation  of  the  vertical  height  at 
which  condensation  will  occur,  which  is  needed  to  determine  the  correct 
radiation  balance.  Thus  fog/cloud  prediction  is  a  natural  consequence  of 
the  model . 


3.  ATMOSPHERIC  BOUNDARY  LAYER  MODEL  CONSIDERATIONS 

The  marine  atmospheric  boundary  layer  (MABL)  is  coupled  both  to  the  overly¬ 
ing  atmosphere  and  to  the  ocean  below.  Thus,  if  a  model  is  to  be  successful , 
it  must  include  large  scale  atmospheric  processes,  such  as  subsidence  and  ad- 
vection,  and  oceanic  processes  which  control  the  sea  surface  temperature. 

To  produce  a  tactically  useful  MABL  model  one  must  consider  the  opera¬ 
tional  environment  in  which  it  will  be  applied  in  addition  to  meteorological 
factors.  Specifically, 

a.  Measurement  of  local  meteorological  parameters: 

Radiosondes  are  normally  released  only  once  or  twice  a  day.  This  in¬ 
formation  is  used  to  determine  atmospheric  conditions  for  the  day.  If  condi¬ 
tions  are  changing  appreciably  it  might  be  necessary  to  launch  another  radio¬ 
sonde  balloon  to  provide  updated  input  information  to  the  model. 

b.  Determination  of  large  scale  parameters: 

In  the  absence  of  data  from  large  scale  models,  a  single  ship  or  sta¬ 
tion  must  be  able  to  calculate  synoptic  scale  forcing,  e.g.  subsidence,  from 
local  measurements.  Several  methods  are  available  to  determine  vertical 
motion  from  radiosonde  soundings,  either  from  one  station  or  multiple  stations 
in  a  region. 

c.  Numerical  prediction  of  local  meteorological  parameters: 

Decisions  must  be  made  in  formulating  a  model  regarding  which  para¬ 
meters  will  be  predicted  by  the  model  and  which  ones  should  be  input  as  fore¬ 
casts  to  the  model  initially.  For  example,  surface  layer  winds  might  come 
from  another  source  and  be  input  to  the  model  when  it  is  initialized,  but  the 
mixed  layer  temperature  would  be  predicted  by  the  model  itself. 
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d.  Possible  alternative  data  sources: 


1)  Ocean  models:  The  atmospheric  and  oceanic  boundary  layers  are 
linked  through  the  sea  surface  temperature  (SST)  which  affects  the  surface 
fluxes  of  heat  and  moisture.  An  alternative  to  measuring  the  SST  would  be  to 
predict  it  with  an  oceanic  boundary  layer  model  coupled  to  the  MABL  model. 

2)  Climatology:  In  the  absence  of  other  information,  climatological 
inputs  could  be  used  for  the  model  to  provide  the  most  probable  boundary  layer 
properties.  Climatologies  are  very  incomplete  at  this  time,  but  should  be 
more  useful  in  the  future. 

3)  Satellites:  Satellite  imagery  might  provide  large  scale  informa¬ 
tion  such  as  SST  patterns.  Much  work  is  being  done  on  extracting  information 
from  satellite  images,  and  their  value,  especially  in  data  sparse  regions,  is 
improving  with  time. 

e.  Output  of  information  in  the  proper  format: 

If  the  output  from  a  MABL  model  is  to  be  used  in  tactical  assessment 
programs  it  must  be  matched  to  the  format  required  by  them.  The  specific  mea¬ 
sure  of  each  quantity  (e.g.  specific  hunidity  or  relative  humidity  for  mois¬ 
ture)  and  the  units  required  by  the  tactical  assessment  program  need  to  be 
determined. 

The  flow  diagram  in  Figure  1  is  for  a  temporal  assessment  of  the  environ¬ 
ment.  (A  spatial  model  would  be  more  complicated  than  Figure  1  shows  since  it 
requires  a  Lagrangian  approach,  special  inputs  of  the  horizontal  variation 
parameter,  and  a  self  consistency  closure  scheme.)  Aspects  of  this  scheme  war¬ 
ranting  further  description  are: 

a.  The  atmospheric  and  oceanic  models  are  considered  separate  entities; 
they  are  linked  through  the  SST  which  is  predicted  in  the  oceanic  model  using 
output  from  the  atmospheric  model. 
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b.  The  propagation  models  use  meteorol ogical  parameters  as  inputs  and 
calculate  propagation  behavior.  Most  of  these  models  do  not  calculate  system 
performance  directly  but  parameters  related  to  the  performance.  The  develop¬ 
ment  of  tactical  aides  that  can  be  used  directly  in  the  field  is  underway. 
Description  of  these  efforts  is  beyond  the  scope  of  this  report. 

c.  The  dynamics  of  the  boundary  layer  are  very  sensitive  to  small  changes 
in  the  meteorological  parameters.  A  sensitivity  analysis  is  needed  to  deter¬ 
mine  the  accuracy  with  which  these  parameters  should  be  specified.  This  will 
impact  on  measurement  systems  and  techniques  and  on  the  design  of  synoptic 


4.  MODELING  THE  BOUNDARY  LAYER 

From  a  local  assessment  perspective,  consider  an  idealization  of  the  ocean¬ 
ic-atmospheric  system.  The  air-sea  interface  is  bordered  by  oceanic  and  atmo¬ 
spheric  turbulent  mixed  layers  (boundary  layers)  which  are  effectively  insu¬ 
lated  from  the  bulk  ocean  and  free  atmospheric  regions.  The  primary  sources 
of  the  turbulence  within  the  atmospheric  (oceanic)  layers  are  the  velocity 
(current)  and  buoyancy  (density)  gradients  near  the  interface.  Even  under  con¬ 
ditions  where  the  water  is  slightly  cooler  than  the  air,  buoyancy  forced  veloc¬ 
ity  fluctuations  within  the  atmospheric  layer  can  be  quite  large  and  mix  the 
entire  MABL  from  the  surface  to  the  inversion.  The  large  vertical  mixing 
yields  constant  (well-mixed)  wind,  temperature  and  humidity  profiles  above  the 
surface.  At  the  top  of  the  atmospheric  mixed  layer  there  is  a  thin  transition 
region  (inversion),  above  it  the  free  atmosphere,  which  is  insulated  from  sur¬ 
face  influence  by  the  inversion.  Idealizations  of  profiles  of  measures  of 
moisture  and  temperature  for  this  situation  are  shown  in  Figure  6. 

The  well-mixed  nature  enables  prediction  of  MABL  evolution  to  be  based 
solely  on  fluxes  at  the  two  boundaries  of  the  layer  (inversion  and  sea  sur¬ 
face)  and  on  the  large  scale  vertical  velocity  and  advection.  This  is  the 
basis  of  recent  model  formulations  by  Deardorff  (1976)  and  Stage  and  Businger 
(1981)  among  others.  Fluxes  at  both  boundaries  are  due  to  buoyant  and  mech¬ 
anically  generated  turbulence.  The  linear  height  variations  of  the  fluxes  of 
the  well  mixed  properties  allow  one  to  relate  fluxes  at  the  inversion  to  the 
more  readily  determined  surface  layer  fluxes  and  general  cloud  features. 

Approaches  exist  for  estimating  synoptic  scale  forcing  from  single  station 
★ 

measurements  but  further  efforts  are  required  to  achieve  the  accuracies 
required  in  MABL  predictions. 

will  be  discussed  in  a  forthcoming  report  on  single  station 
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Physical  processes  in  the  atmosphere  that  control  the  properties  of  the 
MABL  can  be  grouped  as  follows: 

Synoptic  Scale  Forcing 

Synoptic  scale  mean  vertical  motion  and  advection  are  external  factors 
influencing  boundary  layer  evolution.  Advection  brings  new  air  into  the 
locality  being  evaluated  so  that  the  evolution  of  the  boundary  layer  at  a 
point  does  not  depend  solely  on  local  features.  Subsidence  brings  warm  dry 
air  from  above,  and  will  decrease  the  boundary  layer  height.  These  processes 
vary  over  time  scales  on  the  order  of  a  day,  changing  as  the  synoptic  situa¬ 
tion  changes.  A  frontal  passage  could  lead  to  rapid  changes  in  the  synoptic 
forcing. 

Surface  Layer  Fluxes 

The  surface  layer  is  a  region  tens  of  meters  deep  immediately  above  the 
sea  surface.  In  this  layer  wind,  temperature  and  water  vapor  content  change 
from  surface  to  well-mixed  values.  Turbulent  mixing  in  this  layer  is  gene¬ 
rated  by  wind  shear  and  buoyancy.  Due  to  the  gradients,  momentum,  heat  and 
water  vapor  are  transported  upward  into  the  mixed  layer.  The  surface  layer 
thereby  acts  as  a  source  region  for  the  well -mixed  layer. 

Entrainment  and  the  Inversion 

The  inversion  is  a  region  of  wind  shear  and  temperature  and  moisture  gra¬ 
dients.  Thus,  it  acts  in  much  the  same  manner  as  the  surface  layer,  being  a 
region  of  important  fluxes  for  the  boundary  layer.  Turbulence  erodes  the  in¬ 
version,  mixing  the  warm  dry  air  into  the  cool  moist  marine  air.  This  flux 
process  is  called  entrainment.  Entrainment  will  warm  and  dry  the  marine  air 
and  will  cause  an  increase  in  the  boundary  layer  depth.  It  will  also  increase 
the  boundary  layer  wind  speed,  if  there  are  higher  winds  aloft. 
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Radiation 


Several  radiative  heating  and  cooling  processes  are  operative: 

a.  heating  of  the  air  by  solar  radiation, 

b.  heating  of  the  ocean  by  solar  radiation, 

c.  cooling  of  the  sea  surface  by  long  wave  radiation, 

d.  cooling  of  the  cloud  top  by  long  wave  radiation, 

e.  heating  of  the  cloud  region  by  solar  radiation, 

f.  exchange  of  long  wave  radiative  energy  between  the  sea  surface  and 
the  cloud  bottom. 

Condensation 

Cloud  formation  will  take  place  if  the  lifting  condensation  level  is 
within  the  mixing  layer,  that  is,  if  it  is  below  the  inversion.  This  is  an 
important  modeling  consideration  since  the  radiation  balance  depends  so 
strongly  on  the  presence  or  absence  of  clouds.  The  model  predicts  the  inver¬ 
sion  height  and  the  lifting  condensation  level  so  that  a  fog/ stratus  forecast 
is  a  natural  consequence.  This  allows  the  predictive  model  to  branch  between 
clear  sky  and  cloud  topped  cases. 

The  computational  flow  of  the  NPS  model  which  incorporates  the  above 
physical  processes  is  illustrated  in  Fig.  7.  The  initialization  is  based 
on  radiosonde  information,  identical  to  that  required  by  IREPS,  at  the 
beginning  of  the  forecast  period  and  the  specification  must  be  made  for  the 
synoptic  scale  subsidence,  surface  layer  wind  and  sea  surface  temperature 
during  the  forecast  period.  Gleason  (1982)  has  described  and  evaluated  our 
procedures  for  estimating  subsidence  from  single  station  data.  Reformul ation 
efforts  are  now  being  conducted  which  should  allow  the  prediction  of  the 
surface  layer  wind,  given  a  geostrophic  wind  forecast,  and  the  prediction 
of  sea  surface  temperature  based  on  forming  an  initial  upper  ocean  layer. 
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Figure  7 


Schematic  of  input,  prescription  and  computing  steps  in  MABL 
prediction. 


The  predicted  gradients  are  the  mixed  layer  depth,  temperature  and 
humidity  and  the  jump  strength  of  the  matter  at  the  inversion.  The  predicted 
temperature  and  humidity  are  used  to  compute  the  surface  fluxes  and  the  cloud 
base,  if  any.  The  surface  fluxes  and  cloud  conditions  determine  the  entrain¬ 
ment  calculations.  These  differ  for  clear  and  cloudy  skies  and  also  for 
stable  and  unstable  stratification  in  the  surface.  Parameters  deferring  EM/EO 
propagation  properties  of  the  boundary  layer  are  calculated  from  the  values 
and  distribution  of  predicted  and  specified  quantities. 

A  30  minute  time  step  is  used  in  the  prediction  calculation  and  it  was 
based  on  representative  mixing  times  for  the  MABL  (Schacher  et  al ,  1982). 
Running  time  for  a  30  hour  forecast  on  a  HP  9845  computer  is  10-12  minutes. 
Different  running  times  arise  for  clear  and  cloudy  sky  periods  because  of  the 
very  extensive  radiation  calculation  in  the  latter. 
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5.  MODEL  VALIDATION 

Validation  of  the  NPS  integrated  model  approach  is  being  done  for  many 
scenarios  with  data  sets  gathered  in  the  Pacific,  Atlantic  and  Indian  Oceans. 
The  results  show  that  the  technique  is  indeed  a  viable  one,  but  these  results 
must  be  considered  preliminary  in  nature.  The  data  were  obtained  on  shipboard 
and  would  normally  be  available  to  the  meteorology  office  on  a  Navy  ship:  wind 
speed,  air  temperature,  sea  surface  temperature,  hunidity  and  radiosonde 
information.  Acoustic  sounder  results  were  also  available  in  some  of  the 
examined  cases,  but  are  not  needed  to  run  the  model, 
a.  Data  Set 

The  data  to  be  presented  here  were  obtained  during  the  Cooperative  Ex¬ 
periment  on  West  Coast  Oceanography  and  Meteorology  (CEWCOM-78)  conducted  west 
of  San  Nicolas  Island,  CA  during  May  of  1978.  Observations  of  both  oceanic 
and  atmospheric  mixed  layers  were  made  from  the  R/V  ACANIA  and  radiosonde 
observations  were  also  taken  at  surrounding  shore  stations.  The  data  to  be 
shown  are  from  a  48- hour  period,  5/19/0500  to  5/21/0500  PST,  when  the  R/V 
ACANIA  was  cruising  slowly  (2-3  knots)  into  the  wind,  returning  to  a  reference 
point  approximately  every  12  hours.  The  general  location  of  the  R/V  ACANIA 
during  the  5/19  to  5/22  period  and  locations  of  surrounding  shoreline  radio¬ 
sonde  sites  appear  in  Figure  8. 

The  period  was  one  of  steady  onshore  flow  caused  by  the  combined 
effects  of  intensification  of  the  Eastern  Pacific  High  and  the  persistence  of 
the  Mexican  thermal  low.  The  only  apparent  change  in  synoptic  scale  forcing 
during  the  period  was  an  increase  in  the  offshore  pressure  gradient.  Such  a 
synoptic  change  would  have  been  predicted  by  the  geophysics  officer  knowing 
that  the  high  was  intensifying.  Satellite  imagery  showed  increasing  uniform 
stratus  coverage  (thin  to  heavy)  during  the  period  with  a  cellular  (broken) 
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coverage  occurring  late  on  the  21st.  The  increase  in  cloud  coverage  is  pre¬ 
dicted  by  the  model.  We  believe  MABL  evolutions  were  primarily  determined  by 
subsidence,  surface  fluxes,  entrainment  at  the  inversion,  and  cloud  influenced 
radiative  transfer. 

b.  Observed  Changes 

Observations  during  the  19-21  May  period  are  presented  in  Figure  9. 
The  half  hour  average  values  of  wind,  temperature  and  hunidity  measured  aboard 
the  R/V  ACANIA  are  shown  in  9c  and  d.  The  acoustic  sounder  trace  from  the  R/V 
ACANIA  is  given  in  9b.  Atmospheric  soundings  are  shown  in  9a  and  are  compo¬ 
sites  from  the  ship  and  land  stations.  Both  slow  and  rapid  changes  in  boundary 
layer  properties  which  occurred  were  tactically  significant,  [as  will  be 
demonstrated]. 

The  coupled  variations  of  boundary  layer  parameters  are  apparent  by 
examining  the  evolution  on  5/20.  During  that  day:  1)  the  MABL  depth  in¬ 
creased  from  300  m  to  500  m,  2)  the  wind  speed  increased  from  5  to  10  m/sec, 
3)  the  air  temperature  decreased  from  15  to  13°C,  4)  the  sea  surface  tempera¬ 
ture  increased  from  11  to  13°C,  then  went  down  to  12°C  in  the  last  five  hours. 
Changes  in  the  parameters  indicate  that  the  depth  of  the  boundary  layer  was 
influenced  by  entrainment.  Air  entrained  from  above  would  be  warmer,  drier, 
and  have  increased  momentum  due  to  the  higher  wind  speed  above  the  inversion. 
This  would  lead  to  an  increase  in  temperature,  a  decrease  in  relative 
humidity,  and  an  increase  in  wind  speed  as  is  shown  in  Figure  9. 

Possible  coupling  between  the  atmosphere  and  ocean  is  seen  in  the  re¬ 
sponse  of  the  SST  to  the  increased  wind  speed.  The  wind  caused  increased 
mixing  in  the  upper  layer  of  the  ocean,  lowering  the  surface  temperature  by 
mixing  in  cooler  water  from  below.  The  ocean  mixed  layer  deepening  was  evi¬ 
dent  in  bathythermography  records  obtained  from  the  R/V  ACANIA.  The  ocean 
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mixed  layer  deepened  from  10  meters  on  20  May  to  50  meters  on  21  May.  Late  on 
the  20th  the  air  responded  to  the  decreasing  SST,  decreasing  1°C  in  three 
hours. 

An  example  of  the  the  tactically  relevant  changes  in  the  boundary 
layer  is  shown  in  Figure  10  with  respect  to  the  modified  refractivity,  M,  for 
four  different  times.  Changes  of  M  profiles  over  the  three  day  period  of  the 
experiment  were  significant.  The  vertical  dotted  lines  delineate  the  duct 
regions.  At  5/19/1700  a  surface  based  duct  extended  to  300  m.  As  the 
boundary  layer  evolved  the  duct  became  elevated,  and  at  5/21/0500  had  a  base 
at  400  m  and  a  top  at  700  m. 
c.  Model  Predictions 

The  NPS  boundary  layer  model  has  outputs  in  two  basic  formats:  pre¬ 
dicted  time  histories  and  predicted  profiles  of  significant  parameters.  Two 
24-hour  model  forecasts  for  the  validation  period  are  shown  in  Figures  11a  and 
12.  Overlays  of  measured  and  predicted  profiles  of  e,  q,  and  M  are  shown  in 
Figure  11.  The  time  histories  in  Figure  11  also  show  measured  values  of  the 
various  quantities,  and  are  computed  in  half-hour  time  steps.  The  model  was 
re-initialized  at  the  end  of  the  first  24-hour  period,  conforming  to  the  pro¬ 
cedure  that  would  be  used  during  operational  conditions.  Overlays  of  measured 
and  predicted  values  of  z- ,  zlcl,  and  zQ  (evaporation  duct),  and  surface  layer 

O 

values  of  Cn  and  e  are  shown  in  Figure  12. 

The  most  significant  conclusions  that  can  be  drawn  from  comparisons  of 
model  predictions  and  observations  are  as  follows: 

1)  (Figure  11a)  The  mixed  layer  values  of  potential  tem¬ 

perature  (e)  and  water  vapor  mixing  ratio  (q)  were  predicted  quite  well. 
Values  of  potential  temperature,  determined  by  averaging  data  from  several 
stations  in  the  area,  were  predicted  to  within  0.3°C  in  the  first  24  hour 


■M  E 

C  O 

QJ  S- 

+J  s_ 
o  i-  a> 
o.  oj  +j  -o 
-t->  0)  <u 

'mu  E  c 
«J  E  •«- 

-  O  m 

O  r-  -t-> 

CO  o 

f-~ - -  O) 

i  o  • — '  <n 


i  to  a>  -m 
-o  a» 

Ur-  tiv 

«t-  £  a 

-r-  0-  QJ 

E  O  O  U 

S_  r—  X 

O  CL  QJ 

C  QJTJ  fl 

<T3  -•->  - - - 

QJ  -i-  (O 

<_)  CO  «T> 

o  o  >> 

a-Pr- 
TJ  E  •<-  >— 
C  O  -O 

<o  o 

O  QJ  3  • 

•i-  i_  x:  xj 
i-  3  aj « 
Q)  4->  QJ  Of 
sz  id  >  cl: 
a.  i-  -r  to  ■ 
CO  QJ  +->  I 

gaiotj 
E  ■ —  C 
+J  Ql  QJ  'f  . 
C  -M  *-  * 


line  shows  the  boundary  of  the  radar  duct. 


period  and  1.4°C  in  the  second  period.  During  the  first  24  hours  the  mixing 
ratio  was  predicted  to  within  0.9  g/kg;  it  was  predicted  to  within  0.7  g/kg  in 
the  second  24  hours.  These  are  the  basic  meteorological  parameters  on  which 
the  tactical  parameters  depend. 

2)  (Figures  lib,  12a|  The  height  of  the  boundary  layer  and  hence  the 
height  of  the  duct  were  predicted  to  within  14  m  after  12  hours  in  the  first 
period  and  within  58  m  after  12  hours  in  the  second  period.  After  24  hours 
the  errors  in  the  first  and  last  period  were  68  and  75  m  respectively.  The 
duct  was  accurately  predicted  to  rise  from  a  surface  base  to  an  elevated  one 
between  0500  and  1700  on  20  May. 

3)  (Figure  12a)  The  model  predicted  the  lifting  condensation  level 
to  be  below  the  top  of  the  boundary  layer,  thus  predicting  clouds  but  not  fog 
for  the  full  period  as  was  observed. 

4)  (Figures  12c  and  d )  The  model  captures  the  surface  layer  values 
in  Cn  and  e  quite  well,  in  addition  to  correctly  predicting  the  trends  in 
both  24  hour  periods.  The  agreement  arises  in  part  because  both  of  these  are 
quite  dependent  on  quantities  which  were  prescribed  with  the  observed  values. 
Cn  depends  on  the  predicted  air  temperature  and  the  prescribed  sea  surface 
temperature.  8  depends  on  the  predicted  relative  humidity  and  the  prescribed 
surface  wind. 
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Observed  and  predicted  MABL  profiles  for  the  5/19/0500  to 
5/21/0500  CtWCJCM-78  period  including  (a)  Q  (pjn  "*)  and  e 
(°C) ,  (b)  M,  (c)  Cn^  (nT2^) ,  and  (d)  total  extinction  co¬ 
efficient,  6  (km-1).  Solid  lines  correspond  to  observations 
and  dashed  lines  to  model  predicted  values.  Ihe  observed 
values  were  calculated  from  the  meteorological  parameters 
measured  at  the  times  indicated. 


Figure  12.  Observed  and  predicted  surface  layer  values  tor  the  VI'.  '05UU 
to  5/21 /050U  C1MXM-/8  period:  (a)  inversion  height  (Zp 
(solid  line)  and  lifting  condensation  level  (Jx.i.)  (dotted 
line),  'lhe  A's  represent  inversion  heignt  determined  from 
radiosonde  soundings  taken  on  the  k/V  ACANIA,  the  X's  are 
inversion  heights  from  ccmposite  soundings.  'lhe  0’s  repre¬ 
sent  corresponding  observed  LOL’s.  (b)  evaporation  duct 
depth,  Zg  (m)  ,  (c)  (m  ^^) ,  and  (d)  total  (aerosol  and 
water  vapor)  scattering  extinction  coefficient,  e  (km-1). 

In  (b)-(d)  X's  are  the  observed  and  solid  lines  are  the  model 
predicted  values. 
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6.  CONCLUSIONS 

We  have  presented  a  description  of  features  of  the  marine  atmospheric 
boundary  layer  (MABL)  which  influence  the  performance  of  systems  dependent  on 
electromagnetic  and  optical  wave  propagation.  We  believe  the  effects  of  atmo¬ 
spheric  properties  on  such  systems  can  be  estimated  given  accurate  measure¬ 
ments.  However,  the  properties  change  in  both  time  and  space  so  predictive 
capabilities  are  necessary.  We  further  believe  that  predictions  should  be 
based,  as  much  as  possible,  on  local  observations,  on  synoptic  scale  forcing 
discernable  by  an  aircraft  carrier  based  geophysics  officer  and  on  calcula¬ 
tions  which  can  be  made  with  a  microcomputer  such  as  those  planned  for  the 
Tactical  Environment  Support  System  (TESS)  used  for  the  existing  IRE  PS  code 
for  electromagnetic  wave  propagation. 

We  have  presented  a  simplified  slab  model  for  predicting  relevant  changes 
of  the  MABL  due  to  wind  and  buoyancy  forced  mixing  and  cloud  radiative 
effects.  It  requires  as  input,  temperature  and  humidity  profiles  from  an 
initial  radiosonde  and  prescription  of  sea  surface  temperatures,  surface  winds 
and  subsidence  for  the  forecast  period.  It  predicts  the  temperature  and 
humidity  below  the  inversion,  the  inversion  height,  the  magnitudes  of  the 
humidity  and  temperature  jumps  at  the  inversion  and  the  lifting  condensation 
l'vel.  From  the  above  predicted  and  the  prescribed  meteorological  properties, 
one  can  calculate  the  dimensions  of  elevated  and  surface  based  radar  ducts, 
the  intensities  of  optical  turbulence  and  the  magnitude  of  optical  extinction 
due  to  molecular  absorption  and  aerosol  scattering.  The  latter  depend  on  the 
empirical  expressions  relating  aerosol  size  distributions  to  wind  speed  and 
relative  humidity. 

The  NPS  physical  model,  the  initialization  and  the  synoptic  scale  pre¬ 
scription  were  evaluated  from  a  48  hour  period  when  a  cloud  topped  MABL  was 


-40- 


deepening  from  250  to  750  meters.  The  conditions  occurring  were  well  docu¬ 
mented  on  the  basis  of  shipboard  and  surrounding  island  and  shoreline  observa¬ 
tions.  The  most  significant  changes  in  conditions  in  the  tactical  environment 
were  1)  an  initial  surface  based  radar  duct  became  an  elevated  duct  at  end  of 
the  first  24  hours  and  remained  elevated  throughout  the  second  24  hour  period; 
2)  an  overcast  to  broken  stratus  layer  persisted  throughout  the  48  hour  period 
which  did  not  lower  to  cause  a  fog  episode  and  3)  the  surface  layer  wind 
increased  steadily  from  5  to  10  m  s”*  during  the  first  24  hour  period. 

The  model  was  run  on  a  microcomputer ,  which  is  aboard  all  carriers,  and 
predicted  the  first  two  tactical  descriptions  quite  well.  The  surface  based 
radar  duct  was  predicted  to  become  elevated  by  the  end  of  the  first  24  hours 
and  remain  elevated  throughout  the  second  24  hours.  The  predicted  and  observed 
dimensions  agreed  reasonably  well  in  view  of  uncertainties  and  averaging  in 
the  observed  radiosonde  data.  This  agreement  in  duct  dimensions  required 
accuracy  in  the  prediction  of  the  inversion  height,  jumps  at  the  inversion  and 
mixed  layer  temperature  and  humidities.  The  predicted  lifting  condensation 
level  was  always  well  below  the  predicted  inversion  height  implying  an  over¬ 
cast  stratus  layer  with  base  ranging  from  100  m  during  the  first  24  hour  per¬ 
iod  and  lifting  to  200  m  during  the  second  24  hour  period.  The  top  of  the 
layer  corresponded  to  the  predicted  inversion  height.  Observed  lifting  conden¬ 
sation  levels  were  always  higher  (factors  of  2)  than  the  predicted.  The  agree¬ 
ment  on  stratus  cloud  coverage,  even  though  the  bases  were  too  low,  required 
reasonably  accurate  predictions  of  the  mixed  layer  temperature,  humidity  and 
depth  (inversion  height). 

The  wind  speed  increase,  the  third  tactical  description,  was  not  a  pre¬ 
dicted  property  in  the  current  model  but  was  probably  influenced  by  the 
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deepening  of  the  layer.  The  model  is  being  reformulated  to  include  this 
prediction  as  well  as  to  couple  the  model  to  an  ocean  mixed  layer  model. 
These  reformulations  will  be  important  in  achieving  bonafide  predictions  of 
the  other  surface  layer  parameters:  optical  turbulence  (Cn  ),  extinction  ( b ) 
and  evaporation  duct  height  (Zg). 

We  believe  the  demonstrated  reasonable  success  from  present  modeling 
efforts  in  the  prediction  of  the  tactically  significant  features  of  the  MABL 
is  encouraging.  This  encouragement  exists  because  it  requires  a  computer  and 
input  which  are  available  to  a  geophysics  officer  on  a  ship  or  at  a  shore 
facility.  Several  required  areas  of  improvements  and  reformulations  have  been 
identified.  They  are: 

1)  improve  single  station  estimates  of  synoptic  scale  subsidence  and 
advection; 

2)  include  surface  layer  wind  prediction  based  on  synoptic  ( geostrophic) 
wind; 

3)  couple  MABL  to  ocean  mixed  layer  model  so  that  coincident  evolutions 
are  predicted  and  feedback  occurs  through  sea  surface  temperature 
changes. 
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